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ABSTRAGT 

We perform 3D relativistic ideal MHD simulations to study the collisions between high-a (Poynting- 
flux-dominated) blobs which contain both poloidal and toroidal magnetic field components. This is 
meant to mimic the interactions inside a highly variable Poynting-flux-dominated jet. We discover a 
significant electromagnetic field (EMF) energy dissipation with an Alfvenic rate with the efficiency 
around 35%. Detailed analyses show that this dissipation is mostly facilitated by the collision-induced 
magnetic reconnection. Additional resolution and parameter studies show a robust result that the 
relative EMF energy dissipation efficiency is nearly independent of the numerical resolution or most 
physical parameters in the relevant parameter range. The reconnection outflows in our simulation can 
potentially form the multi-orientation relativistic mini-jets as needed for several analytical models. We 
also find a linear relationship between the a values before and after the major EMF energy dissipation 
process. Our results give support to the proposed astrophysical models that invoke significant magnetic 
energy dissipation in Poynting-flux-dominated jets, such as the internal collision-induced magnetic 
reconnection and turbulence (ICMART) model for GRBs, and reconnection triggered mini-jets model 
for AGNs. 


1. INTRODUCTION 

The energy composition in the jet/outflow of astro- 
physical systems is an important and fundamental ques¬ 
tion, since it leads to significant differences for the sub¬ 
sequent energy dissipation process, particle acceleration 
mechanism, radiation spectrum and light curve, polariza¬ 
tion behavior, neutrino emission luminosity, and so on. 
Generally speaking, jets can be separated into two types 
depending on their energy composition: Poynting-flux- 
dominated (PFD) (cr >> 1) and matter-flux-dominated 
(MFD) (cr << 1), where a is the magnetization parame¬ 
ter defined as the ratio between the electromagnetic field 
(EMF) energy flux to the plasma matter energy flux. 

Many independent observations from Gamma-Ray 
Bursts (GRBs), Active Galactic Nuclei (AGNs), micro¬ 
quasars, and Crab nebula give strong hints of the PFD 
outflows at least for some events. Several theoretical 
models have been proposed within the framework of PFD 
jets/outflows to interpret the observations. 

In the field of GRBs, evidence of PFD jets has been 
collected independently in several directions. First, a 
prominent thermal emission compone nt as expected in 
the fi reball-internal-shock model (e.g. Meszaros & Rees 
2000) has been seen o nly in a small fraction of GRBs 
e.g. GRB 090902B, Ryde et al. [2010 Zhang et al 


2011). The majority of GRBs either show no evidence 
of a thermal compon ent or a weak, sub-dominant ther 


mal component (e.g. [Abdo et ah 2009 Guiriec et al 


20111 lAxelsson et al. f2012|. I'hese GRBs require that 
the GRB central engine is highl y magnetized, and jet is 
still PFD at the emission site (Zhang and Pe’er 2009 
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Gao & Zhang 2015). Next, strong linear polarization 
was discovered during the prompt gamma-ray eni ission 
phase for some GRBs ( Yonetoku et al.pOll 2012), and 
during the reverse-shock-dominated early optical after¬ 
glow emission phase for some others (Steele et al.|2009 
Mundell et al.j20T3 ), which hint at the existence of glob¬ 


ally ordered magnetic fields in the jet. Furthermore 
strong PeV neutrino emission as predicted by the MFD 


models has not been observed from GRBs so far (Ab- 


basi et al. 2012), which is consistent with the expecta- 
tion of the PL'D models (jZhang & Kumar|2013|). Finally, 
the MFD internal shock (IS) model for GRBs also suf- 


ciency (Panaitescu et al.|1999l 

Kumar|1999), electron fast 

cooling ( 

Ghisellim et al.ipOOO 

1. the electron number ex- 

cess (By. 

kov ^ Meszaros||1996 

jDaigne and Mochkovitch 

1998 Shen and Zhang 2009 

), and inconsistency with 

some empirical (Amati/Yonetoku) relations (Zhang and 


Meszaro^|2002 Liang et al.||2010| ). Zhang Van ( |201l[ ) 


proposed a novel PL'D outflow model named as “the 
Internal-Collision-induced MAgnetic Reconnection and 
Turbulence (ICMART)”, which can potentially keep the 
merits of the IS model but alleviate the criticisms faced 
by the IS model mentioned above. The main idea of 
the ICMART model is that the GRB jets are Poynting- 
flux dominated. The Poynting flux is catastrophically 
discharged at a relatively large distance (e.g. 10^^ cm) 
from the central engine through collision-induced mag¬ 
netic reconnection. The magnetic energy is converted to 
particle energy and radiation efficiently, leading to a very 
high radiation efficiency as demanded by the G RB data 
(jPanaitescu & Kumar||2002| jZhang et al.||2007|). A PFD 
jet has less leptons than the ML'D model so that the elec¬ 
tron excess problem is avoided. A large emission radius 
favors a moderately fast cooling, which can account for 
the right low-energy spectral index observed in GRBs 
(Uhm & Zhang 2014). It also gives a natural explana- 
tion of the seconds-duration of “slow variability compo- 
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nent”observed in GRBs (Gao et al. 2012). The rapid 
“fast variability component” can be interpreted within 
this scenario as r nini-jets due to locally Lorentz boosted 
regions (see also Lyutikov & Blandford 2003 Narayan 
& Kumar||200LM. It is speculated that turbulent recon¬ 
nection in a moderately high-a flow can give rise to rel¬ 
ativistic motion of mini-jets within the bulk relativistic 
motion of the jets. 

For AGNs, observations show fast variable TeV flares 
of two blazars (Mrk 501 and P KS 2155-304) (Aharonian 
et al.||2007t |Albert et al.||2007|). The light-crossing time 
IS even shorter than the event horizon size of the black 
holes, so that emission must come from a small local re¬ 
gion. The derived Lorentz factor in the emission r egion 
should be larger than 50 (Begelman et al. 2008 Mas- 


tichiadis & Moraitis 2008). 'I'his value is much larger 
than the observed Lorentz factor of the bulk motion of 


emission. Although these models show great potential to 
interpret the observations and overcome the criticisms in 
the traditional MFD models, some important ingradients 
of the models are still of a speculative nature. Detailed 
numerical simulations are needed to give a solid footing 
to these models. 

From the morphologic point of view, jets/outflows can 
be categorized into two types: continuous and episodic. 
Theoretically, episodic jets can be formed either from 
a highly variable central engine with variable accretion 
rate; or dis ruption of a continuous jet by screw or kink in¬ 
stabilities ( Li||2QQQ Mizuno et al.||2QQ9 ); or from a MHD 
erruption process s imilar to solar corona l mass ejection 
(|Yuan et al.||2009t |Yuan & Zhang||2012|). Observation- 
ally, episodic jets or knots in jets have been observed in 
many X-ray binaries ([Mirabel fc Rodrignez|1994 [Hjellm- 


the global jet, which is generally smaller than 10 (Giro- 
letti et al.| 20041 |Pi ner fc Edwards]|2004|). To interpret 
these observations 


ing fc Rupe 


naries (iMirabel <!v Kodrignez|P 
1 995| ^p dgL_^ J^^llQ^il 


Doi et al. 


llVlar scher et al.||2002l |Chatterjee et al.J|2009 
01 1[ ). Rapid variabilities ob served in GR'Bs also point 


Gianmos et al.l (|2009p proposed a towards highly episodic jets (Rees and Meszaros 1994 


“jets in a jet” model, which considers that some mini-jets 
are generated by local reconnection outflows in a global 
PFD jet. The mini-jets can give extra Lorentz boosting 
and particle acceleration to generate the observed TeV 
photons around these loc al reconnection regions with fast 
variability. Even though Giannios et al. (2009) did not 
specify the mechanism of magnetic dissipation, observa¬ 
tions of AGN jets reveal bright knots that are consistent 
with internal interactions within the jet. Within the PDF 
jet scenario, IGMART processes similar to what are en¬ 
visaged in GRB jets may also play a role. 

Another related astrophysical phenomenon is y-ray 
flares observed from the Grab nebula. Monte carlo sim¬ 
ulations suggest that the bright y-ray flares and fluctu¬ 
ations in longer time scales can be understood within 
the framework that there are many mini-jets with a wide 
distribution of size and Lorentz factor within the PFD 
outflow of the pulsar. The flares correspond to the e pochs 
when some bright mini-jets beaming towards earth (Yuan 
et al.|2QlT ). The observations suggest that similar physi¬ 
cal processes as those operating in GRBs and AGNs may 
be playing a role in the Grab nebula. 

In ano ther front, recent Partical-In-Gell (PIG) sim- 


ulations (Sironi & Spitkovsky 
show that reconnection under 


Mil |Guo et al-llMi ) 


ligh-a condition can ef 


flciently accelerate thermal particles to form a non- 
thermal power-law population of the particles. This gives 
a good support to the above PFD models from the par¬ 
ticle acceleration point of view. 

The models discussed above for different astrophysi¬ 
cal systems share some common physical processes, such 
as efficient magnetic energy dissipation in the PFD out¬ 
flow/jet, mini-jets generated by the relativistic outflows 
due to local reconnections, particle acceleration in the 
reconnection region, and production of the non-thermal 


^ |Lyutikov & Blandford| (|2003| and jNarayan & Kumar | (|2009| 
proposed that UK-K variability is a consequence ot mini-jets due 
to relativistic outflow from reconnection or relativitic turbulence. 
There is no simple explanat ion to the observed slow variability 
component in these models. |Zhang & Yan| (|2011|) attributed the 
two variability components (slow and tastj as due to central engine 
activity and mini-jets , respectively. Monte Carlo simulations by 
|Zhang &: Zhang|(|2014|) showed that the ICMART model can indeed 
reproduce tiie observed GRB light curves. 


Paczynski and Xu 1994). As a result, studying iriterac 
tions or collisions between magnetic blobs or shells is of 
great interest. 

In this paper, we perform detailed numerical simula¬ 
tions on the global properties of collisions between high- 
( j blobs, as envisaged in the IGMART model of GRBs 
(jZhang & Yan||2011|). In Section p) we give a brief 
introduction of our 3D relativistic MHD code and the 
simulation setup. In Section ^ we present an exam¬ 
ple simulation case to show the Key results, and perform 
a detailed analysis and resolution study. We then ex¬ 
pand our simulations on two-blob collisions in Section 
Q to a large parameter space and discuss how differ¬ 
ent parameters affect the simulations results. In Section 
we show preliminary results for multiple collisions 
among four high-cr blobs. We summarize our results in 
Section ^and discuss the implications of our simulation 
results on some high energy astrophysical systems, such 
as GRBs and AGNs. 


2. NUMERIGAL METHOD AND PROBLEM SETUP 
2.1. Code introduction 

We use a 3D special relativistic MHD (SRMHD) 
code which solves the conservative form of the ideal 
MHD equations using higher-order Godunov-type finite- 
volume methods. This code is a development version of 
the “LA-GOMPASS” MHD code which was first devel¬ 


oped by Li & Li (2003) at Los Alamos National Labora- 
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tory. The equations solved in the code are: 
d{Tp) 


V 


dt 
aC C 2 


(TpV) = 0, 

E X B 


Attc 

E 0 E + B iSi B 


r^h 

) + v-[^v< 


‘V + {p + 




( 1 ) 

^2 


Stt 


An 




■p - rpc^ + 


= 0 , 

B'^+E'^ 

Stt 


( 2 ) 


) + v ■i(r^h-rpc^)v 


+ 7-® X B] = 0 , 

47r 

(3) 

dB 

(4) 

+ cV X P — 0, 
dt 

V 

E =-X P, 

c 

(5) 


where T, p, h, P are the Lorentz factor, rest mass density, 
relativistic enthalpy, and gas pressure, respectively, V, 
E, B are the vectors of fluid velocity, electric field, and 
magnetic field, respectively, and the symbol “(8)” denotes 
tensor product. We also use the ideal gas equation of 
state: p = (7 — l)u^ where 7 and u are the adiabatic 
index and the internal energy density, respectively. 

We use HLL flux with the piecewise paraboli c recon¬ 


struction method to solve the Riemann problem (Colella 


& Woodward 1984), and use the constrained transport 
(CT) method (jbalsara & Spicer||1999t |Guan et al.||2014 ) 
to ensure V • = U. We use the Cartesian coordinates 

(x^y^z) in our simulations. 


2 . 2 . Problem set up 

We envisage that the central engine of GRBs or AGNs 
launch a Poynting-flux-dominated jet/outflow. As dis¬ 
cussed in Sect. 1 , episodic jets are preferred from ob¬ 
servational data. Even if the jet may be overall continu¬ 
ous, it is very likely non-uniform internally and may form 
many knots in the jet, where a much larger amount of 
EMF energy (Pem) is concentrated compared with other 
sparse regions in the jet. We can simplify the knots of the 
jet/outflow as many quasi-isolated magnetic blobs with 
both poloidal and toroidal field components. Due to the 
intrinsic erratic behavior at the central engine, different 
magnetic blobs may have different velocities at the emis¬ 
sion region, so that multiple collisions are very likely to 
happen among different blobs. Due to the ultra rela¬ 
tivistic motion of the jet, the relative velocities between 
different blobs can easily become relativistic. _ 

I n our s imulation domain, we use the model from |Li et] 
ah] (|2006|) to initialize the magnetic field configuration. 
The equations are introduced in the cylindrical coordi¬ 
nates (r, 0, z), and we will transfer them to the Cartesian 
coordinates (x, p, z) in our simulations, from the center 
(r = 0 ) of each blobs, the field is assumed to be axisym- 
metric with the poloidal flux function 4> as 


^{r,z) = Bb.or^exp ( - 


r^ + z^ 


(6) 


and the relationship between 4>(r, 2 :) and the 0 compo¬ 
nent of the vector potential is 4>(r, z) = rA^. Pb,o 
and To are the normalization factor for the magnetic 
strength and characteristic radius of the magnetic blob. 


respectively. One can then calculate the r— and 2 :— com¬ 
ponents of the poloidal field 



The poloidal field is closed and keeps the net global 
poloidal flux as zero. The toroidal field configuration is 
motivated by considering the black hole accretion disk 
system as a “dynamo”, which shears the poloidal flux to 
form the toroidal flux from the rotation. The toroidal 
component of the magnetic field therefore has the form 


B(f) 


ror 


r f + z^ 

Bb,oa— exp- 5 — 

ro V 


(9) 


Here the p arameter a cont rols the toroidal-to-poloidal 
flux ratio. Li et al. (2006) showed that when a ^ 3, 
the two flux components are roughly equal with each 
other. We set a = 3 for our example sim ulation, and 
explore a larger value of a in Section §4.7| We choose 
the comoving center-of-mass frame of the blobs as our 
simulation frame. The direction of velocity is along Z- 
axis with a profile 



^,Z5 ('^ ^ ^ 0)5 

Vb,zexp ^ , (r > ro), 


( 10 ) 


where Tb,z is a constant value which can be either positive 
or negative corresponding to -\-Z or —Z direction of the 
velocity. We also set a uniform gas pressure value (P) 
both inside and outside the blobs. The value of P is 
much smaller than the initial magnetic energy density of 
the blobs. 

For the density profile, we first define a constant initial 
value of the blob magnetization parameter around the 
central region of the blobs: 




Pem 

T^’ 


( 11 ) 


where h = pc^ + jP/il — 1 ) is the specific enthalpy de¬ 
fined in the fluid’s comoving frame, p is the rest mass 
density, P is the gas pressure introduced above, 7 is the 
adiabatic index, T is the bulk Lorentz factor calculated 
by the velocity profile introduced above, and Pem is the 
EMF energy density calculated by Pem = + P^)/ 87 r 

from the magnetic held profile introduced above. The 
density profile is therefore 


P = 


( 


-^em 

r2ab,i 


tP 

7-1 



r < ro 


and p > pbkg), 


Pbkg5 


(r > ro or p < pbkg), 


( 12 ) 

where pbkg is a constant parameter to control the uniform 
background mass density. 

We also introduce two position-control parameters Zd 
and Xg. For a collision between two blobs, the center of 
the two blobs are located at (xi, zi) and (x 2 , ^ 2 ), so 
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= \zi — Z 2 \ is the initial distance between the center 
of the two blobs in Z direction, and Xs = \xi — X 2 I is the 
initial misalignment between the center of the two blobs 
in X direction due to the possible misalignment of the 
blobs. The Y coordinate is the same for both of them. 

In Table we give the normalization relationship be¬ 
tween the code units and the physical units. There are 
only three free parameters, Lq, Bq and c to control the 
normalization of the entire system. Defining different 
physical values of Lq and we can normalize the sim¬ 
ulation system to different environments and problems. 
In Table we also list two sets of example typical val¬ 
ues to show the way of the normalization. In the rest of 
paper, all the parameters are given using code units. We 
keep vq = 1.0 for all the following simulations. In addi¬ 
tion, we use 7 = 5/3 in most of the simulations, since 
most of the regimes are mildly relativistic. This may not 
always be true, especial ly in the regions of reconnection 
outflows, so in Section ^4.8 we also try 7 = 4/3 to test 
the difference. 


3. AN EXAMPLE CASE 

In this section, we show a series of detailed analyses 
based on one example simulation case. We focus on the 
following aspects: the evolution of magnetic energy to 
address the efficiency of magnetic energy dissipation, the 
details of the collision process, the properties of magnetic 
reconnection and outflows, and the numerical resolution 
effects. We reveal significant collision-induced reconnec¬ 
tion events with a remarkable efficiency around 35%, 
which is resolution insensitive. The outflow properties 
of reconnection events indicate the potential capability 
to generate super-Alfvenic relativistic mini-jets. 

3.1. Initial parameters 

The initial parameters of the example run are listed 
in Table [21 We consider two identical blobs with initial 
magnetization parameter crb,i = 8 separated by = 4.4, 
with an X-direction offset 1.0. The two blobs move in 
opposite directions in Z direction with an initial center 
speed Eb,z = 0.3 c. The background pressure and den¬ 
sity are P = 10“^ and pbkg = 10“^, respectively. In 
order to clearly show the initial magnetic field configu¬ 
ration of the blobs, in Eigure we show a ^ = 0 slice 
(cut through the blob centers) of the profiles of several 
parameters: projected field line configuration (panel A), 
a distribution (panel B), B^ (panel C), and By (panel 
D). 

Eor this example run, the 3D box size is chosen as 20^ 
from -10 to +10 in each dimension, which means that 
the position (x,y,z) = ( 0 , 0 , 0 ) corresponds to the center of 
the box. And the resolution is chosen as 1024^. 

3.2. Energy evolution analysis 

Since the initial magnetic configuration is not in com¬ 
plete force balance (between the internal magnetic pres¬ 
sure and the background gas pressure), the blobs would 
quickly expand and evolve into a quasi-steady phase, 
forming a quasi-force balance between the gas pressure 
and magnetic pressure. During this process, a fraction 
of EME energy E^em is converted to thermal and kinetic 
energy due to magnetic field relaxation. So before per¬ 
forming a collision simulation, we first simulate the blob 


evolution of non-collision case to quantify the EME en¬ 
ergy level in the quasi-steady phase. This would serve 
as the reference value to be compared with the collision 
case in which additional EME energy drop is expected 
due to additional magnetic dissipation. 

The upper panel of Eigure shows the evolution of the 
blob electromagnetic energy Eem as a function of time 
(normalized to the initial value Eem,o)- The dashed line 
shows the evolution in the non-collision case. There is 
a significant drop of Eem before t ^ 6 , which is due to 
the magnetic field relaxation during the process of es¬ 
tablishing a force balance between the outward magnetic 
pressure force and the inward gas pressure force. After 
the balance is established, Eem is nearly constant and 
enters a quasi-steady phase, which can be used as the 
reference energy level without collision. 

Next, we simulate the collision case between two high¬ 
er blobs. The initial parameters for these two blobs are 
the same as the non-collision case. The Eem evolution of 
the two blobs with collision is shown as the solid line in 
the upper panel of Eigure The efficiency ( 7 ) of 
energy dissipation due to collision-induced process can 
be calculated by 


T] = 


T'errijnc Bem,c 


Ee 


(13) 


where Eem,c and Eem,nc are the EME energy values for 
the collision and non-collision cases, respectively. The 
efficiency of the example case is shown in the lower panel 
of Eigure where we find that the efficiency is about 
35% near the end of collision process. This efficiency 
is much higher than the collision-induced kinetic energy 
release efficiency in the MED outflows in the internal 


or less (e., 

g. Panaitescu et al.|1999 

Kumar|1999 

Maxham 

& Zhang 

20U9[ Gao IVIeszaros 

2U14|). It is consistent 


& Yan 2011 see more discussion below in ^3.5) 


One important question is what mechanism causes this 
efficient magnetic energy dispassion? Erom the magnetic 
configuration we can see B^ and By have opposite di¬ 
rections around the collision region (see Eigure]^. We 
suggest that most likely the additional E^em dissipation is 
triggered by strong collision-driven reconnection events. 
In order to check our conjecture, in the following, we 
carry out a series of detailed analyses based on our sim¬ 
ulation data. 

The E^em evolution in Eigure can be characterized in 
four stages: ( 1 ) an initial “self adjustment” (steep decay) 
phase before t ^ 10 ; ( 2 ) a “plateau” phase from about 
t ^ 10 to t ^ 38; (3) a “normal decay” phase from about 
t ^ 38 to t ^ 120; and (4) a final quasi-steady phase. We 
analyze these stages in detail below. 

The major collision starts from the later part of the 
“self-adjustment” steep decay phase. The collision com¬ 
presses the magnetic fields to make the energy level 
higher than non-collision case. Panel A of Eigureshows 
a series of representative cuts at t = 4. Erom left to 
right, the four images display the 3D current contour plot 
viewed from Y-axis, the 3D current contour plot viewed 
from X-axis, the 2 D contour cut of the ^-component of 
the outflow velocity (Vy) in the YZ-plane (x=0) corre¬ 
sponding to the current plot, and the 2D contour cut 
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TABLE 1 

The normalization factors between physical units and code units. 


Parameters: 

Length 

Velocity 

Time 

Magnetic field 

Pressure 

Density 

Code units: 

1 

1 

1 

1 

1 

1 

Normalization factors: 

Lq 

c 

Lo/c 

Bo 


Bllc^ 

Typical values 1: 

10^^ cm 

3 X IQio cm/s 

33 s 

10^ G 

10® Ba 

1.1 X 10-1^ g/. 

Typical values 2: 

cm 

3 X IQio cm/s 

333 s 

10 G 

102 Ba 

1.1 X 10-19 g/. 



Panel C: Bx 



Panel D: By 



Fig. 1. — Several manifestations of the initial magnetic field configuration cut in the blob-center plane in the example simulation. Panel 
A: The initial 3D field line profile viewed along the Y direction. The color contour denotes the value of Panel B: The 2D contour 

cut of the initial a profile in the XZ plane (y = 0); Panel C: The 2D contour cut of the ar-component of the initial magnetic field strength 
in the XZ plane (y = 0); Panel D: The 2D contour cut of the y-component of the initial magnetic field strength in the XZ plane (y = 0). 


TABLE 2 

The initial parameters for the example simulation. 



Bb,o o: 

l^b,z| 

P 

Pbkg Zd Xs 

8 

VTr 3 

0.3c 

10-2 

10-1 4 4 0 


of the rest mass density in the FZ-plane (x=0), respec¬ 
tively. From these results we find that a strong current 
layer and a pair of outflows are forming around the con¬ 
tact surface, which are consistent with the features of a 
collision-driven reconnection. 

The second stage is the “plateau” phase from about 
t ^ 10 to t ^ 38. Panel B of Figure shows a se¬ 
ries of representative cuts at t = 18. We can see that 
the current layer around the contact surface becomes 
clearer and more concentrated. The outflows become 
faster (nearly 0.75c) and are also more concentrated 
at the current layer. Besides the four representative 
cuts shown in all panels, for panel B, we also add one 


extra 2D contour cut of the x-component of the out¬ 
flow velocity (14) in the XZ-plane (y=0) corresponding 
to the current plots, which presents another important 
result that the current layer actually generates multi¬ 
orientation outflows in a 3D structure. These results sug¬ 
gest that many mini-jets with relativistic speeds can be 
potentially generated, if multiple collisions are invoked 
in a PFD outflow. Another interesting phenomenon is 
that although the system undergoes a strong reconnec¬ 
tion process which in principle dissipates the EMF energy 
significantly, the global F^em evolution is nearly fiat and 
even shows slight increase during this stage. The main 
reason for this feature is that the initial strong recon¬ 
nection is collision-driven. Besides the strong reconnec¬ 
tion, collision-induced strong compression also exists and 
tends to increase which balances and even slightly 
surpasses F^em dissipation due to reconnection. Th e addi¬ 
tional outflow study in the following Section §3.3[ which 
shows that the outflows become super-Alfvenic at this 
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Fig. 2. — Upper panel: The Poynting flux energy (-Fem) evolu¬ 
tion of the example simulation case. Dashed line denotes the non- 
collision case, which serves as the reference for additional magnetic 
dissipation. Solid line denotes the case of collision between two 
blobs, lower panel: Ratio calculated by (-Fem,nc — -E^em,c)/-E^em,nc 
to show the additional -Fem dissipation efficiency triggered by the 
collision-induced processes. 


stage, also supports the above analysis. 

The next stage is the “normal decay” phase. We choose 
two series of representative cuts at t = 58 (Panel C) and 
t = 94 (Panel D), respectively. The current strength 
and outflow velocity are similar between panels C and 
D, while they are systematically weaker and slower com¬ 
pared with the “plateau” phase (panel B). This means 
that the initial collision-driven effect becomes weaker 
and the reconnection-facilitated dissipation enters a rel¬ 
atively steady phase. In the mean time, compression 
becomes sub-dominant, so that globally L^em dissipates 
with a relatively steady rate, which roughly equals to 
^ ~ center-of-mass frame of the 

blobs (I/n is the length normalization factor introduced 
in Table Q . The additional outflow study in the follow¬ 
ing Section §3.3[ which shows that the outflow velocity 
keeps being around the Alfven velocity at this stage, also 
supports this conclusion. 

Finally, after t ^ 120, the reconnection-dissipation 
gradually becomes weaker, and the system enters the 
quasi-steady phase without obvious F^em dissipation. The 
F^em evolution becomes nearly parallel with the non- 
collision case in Figure 

From these analyses, we conclude that the collision be¬ 
tween two high-cr blobs can indeed trigger strong mag¬ 
netic reconnections and dissipate a significant fraction 
of EMF energy due to the reconnection-facilitated pro¬ 
cesses. 


3.3. Additional outflow study 

Following the above analyses, in this subsection, we 
carry out an additional study on the outflow velocity. 
We compare the local Lorentz factor of the outflow (Font) 
with the critical Lorentz factor Fa calculated from the 
local relativistic Alfven velocity 

Ia = , (14) 

UIirh/JWTT 

and the critical Lorentz factor Fms calculated from the 
maximum possible value of the local relativistic fast 


mode magnetosonic velocity 

= ^Fi + Cf(l-Vi/c2), (15) 

where h' and B' are the specific enthalpy and magnetic 
strength in the local comoving frame of the fluid, and Cg 
is the relativistic sound speed calculated by 


Cs = c^/WJh'- 


(16) 


In order to investigate whether the fluid velocities ex¬ 
ceed the two characteristic velocities, we define 


Ra = 


Bms — 


J- out 

T7’ 

bout 


(17) 

(18) 


Figure shows the selected 2D contour cuts of Rp^. 
The three panels in the upper row correspond to the 
starting time when Rp^ > 1 is reached (t = 4), the time 
when i?A is the largest (t = 18), and the ending time 
for the condition of Ra > 1 (t = 38), respectively. After 
t ^ 38, the Tout starts to become slightly smaller but 
still close to Fa (see the three panels in the lower row of 
Figure]^. These results are consistent with t he e nergy 
evolution analysis presented above in Section ^3.2 The 
duration when Ra > 1 is satisfied is just the “plateau” 
phase of energy evolution, in which strong compression 
exists and drives the outflows to become super-Alfvenic. 
After t ^ 38 the energy evolution enters the “normal 
decay” phase, which corresponds to the phase of rela¬ 
tively steady reconnection-facilitated dissipation without 
strong compression, so that the outflow velocity is close 
to the theoretical Alfvenic velocity. 

Figure show the contour cuts of i?ms- Since Wis is 
the maximum wave propagation speed in a MHD sys¬ 
tem, if i?rns > 1, a local shock in the front of the out¬ 
flow would potentially be generated. The three epochs 
shown in Figure [^correspond to the starting time when 
i?ms > 1 is satisfied {t = 6), the time when i?ms is the 
largest (t = 18), and the ending time for the condition of 
i?ms > 1 to be satisfied (t = 20), respectively. These re¬ 
sults indeed show a period of about 15 time units during 
which i?ms > 1 is satisfied. This duration is shorter than 
the duration when Ra > 1 is satisfied. For this case, 
the largest value of i?ms is about 1.13. Sinc e the Tout de¬ 


pends on numerical resolution (see Section ^3.4 below for 
details) and other physical parameters, it is worthwhile 
to perform a more detailed study for this feature in the 
future. In this study, since Vbut is only slightly larger 
than Vms in a small local region and for a short dura¬ 
tion, we do not resolve an obvious shock feature from 
the numerical data. 


3.4. Resolution study 

We now discuss the effects of numerical resolution on 
our results. Although the ideal MHD code that we use 
does not have explicit resistivity, it still has numerical 
resistivity from the numerical scheme, which depends on 
the resolution of the simulation. This may affect the 
reconnection rate and energy dissipation rate in the sim¬ 
ulations. To address this uncertainty, we perform a reso¬ 
lution test based on the above example case. We keep the 
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Panel A: stage 1, cut at t=4 



Panel B: stage 2, cut at t=18 




l0.135 

I0.I2O 

I0.IO5 


0.090 

0.075 

0.060 

1 0.045 

0.030 

0.015 


Panel C: stage 3, cut at t=58 


Current ' I Current - 



_► 0 I _ 0 

^ ^ ^ h " ^ 



Panel D: stage 4, cut at t=94 




Fig. 3. — The representative cuts of current, velocity and density for the different evolution stages corresponding to Figure Panel A 
corresponds to the initial “self adjustment” phase; Panel B corresponds to the following “plateau” phase; Panel C & D correspond to the 
“normal decay” phase. The last quasi-steady phase has no obvious feature, so we do not draw cuts for that stage. For each panel, the cuts 
from left to right are the 3D current contour plot viewed from T-axis, the 3D current contour plot viewed from A-axis, the 2D contour cut 
of the y-component of the outflow velocity (Vy) in the TZ-plane (x=0) corresponding to the current plot, and the 2D contour cut of the 
rest mass density in the TZ-plane (x=0), respectively. In Panel B, we add an additional 2D contour cut of the ^-component of outflow 
velocity (Vx) in the AZ-plane (y=0) to show the existence of multiple directions of the outflows. 

same box size and the parameters in Table and only to the results with numerical resolution 128^, 256^, 512^ 
change the resolution. Figure shows the results. The and 1024^, respectively. When the resolution decreases, 
magenta, red, green and blue groups of lines correspond we find that the level of F^em evolution is systematically 
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Fig. 4.— The selected 2D contour cuts of for different stages. The three panels in the upper row correspond to the starting time 
when Ra > 1, the time when Ra is the largest, and the ending time for the condition of Ra > 1, respectively. These correspond to the 
plateau stage. The three panels in the lower row correspond to three epochs during the normal decay phase, during which Pout becomes 
relatively steady and close to Fa- 



Fig. 5.— the selected 2D contour cuts from the results of Rms- The three panels correspond to the starting time when Rms > 1 is 
satisfied, the time when Rms is the largest, and the ending time when the condition of Rms > 1 is satisfied, respectively. 


lower and the efficiency also slightly decreases. On the 
other hand, the change of efficiency is only several per¬ 
centage from the highest to the lowest resolutions, which 
means that the L^em dissipation efficiency is insensitive 
to numerical resolution. In addition, the L^em level and 
the efficiency in the final quasi-steady phase also show a 
trend of convergence when the resolution increases. 

Another important result from the resolution study is 
that the maximum outflow velocity increases when the 
resolution increases. Figureshows the contour cuts cor¬ 


responding to the maximum ^-component of the outflow 
velocity (Vy) in the YZ-plane (x=0) for different resolu¬ 
tions. The maximum values of Vy are about 0.45c, 0.55c, 
and 0.75c for resolution of 256^, 512^, and 1024^, respec¬ 
tively. The reason is probably that the higher resolution 
decreases the effective numerical resistivity and decreases 
the aspect ratio between the thickness and the length 
of the reconnection layer, so that the outflow speed is 
forced to reach a higher value in order to balance the 
similar compression forced inflow. This analysis is also 
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Fig. 6.— A numerica l res olution study based on the above ex¬ 
ample case in Section §3.2| The magenta, red, green and blue 
groups of lines correspond to the resolutions of 128^, 256^, 512^ 
and 1024^, respectively. The -Fem dissipation efficiency at the finial 
quasi-steady phase is nearly the same in all cases. 


supported by Figure which shows that with an in¬ 
creasing resolution, the length of the reconnection layer 
is similar, but the thickness becomes thinner. 


3.5. Physical analyses 

In this subsection we carry out some physical analy¬ 
ses to understand the ^ 35% F^em dissipation efficiency 
obtained from our numerical simulations. 

Assuming a complete inelastic collision between two 
high-<j blobs, Zhang & Yan (2011) analytically estimated 
the total efficiency of the collision-induced ^em dissipa¬ 
tion efficiency based on energy and momentum conser¬ 
vation laws. Their Equation (51) can be written as 


T] = 


^m{mi + 7742 ) 


1 + Crb,f (FlTTli + r27742)(l + Crb,i) ’ 


(19) 


where crb,i is the initial a value of the two colliding blobs, 
(Jby is the final a value after the inelastic collision is 
over, Fi, r2, and F^ are the Lorentz factors of the two 
colliding blobs and the merged blob, respectively, and 
774i, 1712 are the masses of the two colliding blobs. In our 
simulations, the two blobs are identical so that ttt-i = m 2 . 
Since we are observing in the merged frame so that F^ = 
1 , Fi = r2 = r, the final expression of the efficiency can 
be reduced tcE] 


1 1 

1 + crb,f r(l + crb,i) 


( 20 ) 


In order to connect this analytical equation with our 
simulated results, we first carry out some analyses to see 
if the condition of complete inelastic collision is satis¬ 
fied. For ideal MHD simulations, fluid elements are at¬ 
tached to the field lines. Tracking the evolution of mag¬ 
netic field configuration is therefore a convenient way to 
study whether collision is inelastic. Figure shows sev¬ 
eral contour cuts of the 3D field line evolution. Initially 
the fields are compressed around t = 6, and then bounce 
back around t = 12. Later strong collision-driven recon¬ 
nections on the contact surface efficiently dissipate the 
compressed magnetic energy and reduce the magnetic 
pressure in the center. This prevents further bouncing 


^ This can be also derived directly by writing energy and mo¬ 
mentum conservations in center-of-mass rest frame. 


back and reorganizes the field configuration to make the 
two blobs merge into one larger blob with a new field 
configuration with a “ 00 ” shape at the final quasi-steady 
stage. This suggests that the two blobs merge to one 
entity after the collision. 

Due to the mis-alignment in X direction of the two 
blobs, the collision would induce rotation during the 
merging process. This would render the collision pro¬ 
cess not completely inelastic. In order to investigate how 
important the rotation effect is, we calculate the ratio be¬ 
tween the rotation energy (E^ot) and the initial kinetic 
energy (F^k,i)- The rotation energy of the two blobs can 
be estimated as E^ot = 2x(|/cj^), where the moment of 
inertia of one blob can be estimated as / = ^mr‘^ 
where the first term denotes the moment of inertia of 
an idea sphere, and the second term denote the displace¬ 
ment from the rotation axis. Since the blobs expand with 
time, the size of the blob and its displacement increase 
with time. We estimate that after the merging process, r 
is about three times of ro = Fq = 1. We therefore derive 
/ ^ ^774 Fq = ^ 774 . For the angular velocity cj, we can 
estimate it from Figure!^ which shows a roughly 7r/4 
angular change within At = 90Lo/c = 90. As a result, 
one can estimate uj ^ 7r/360, so that Erot = 10“^774. The 
ratio between Erot and Fk,i is therefore 


En 


Ek,i 


10 ^774 


10 


-2 


( 21 ) 


So the rotation energy is only a small fraction of the 
initial kinetic energy, which means that the collision is 
very close to completely inelastic collision for this ex¬ 
ample case with Xg = 0.5. While even if F^ot becomes 
a larger fraction of Fk,i when the misalignment Xg in¬ 
creases, it would only reduce the kinetic energy dissi¬ 
pation efficiency, but would have little direct effect on 
the Fem dissipation efficiency that is our primary con- 
cer ■13 Due to the initial high-cr property of the blobs, 
the contribution from the Fk,i dissipation to the total 
dissipation efficiency is only a minor fraction when Eem 
has significant dissipation, as we have found above. 

With the above preparation, we can achieve a phys¬ 
ical understanding of the high efficiency obtained from 
our simulation. Based on Eq.(20), we can derive the ex¬ 


pected efficiency. Erom the initial condition, we derive 
F = 1.05. From simulation results, we can also calcu¬ 
late crb,f. Since dby has a complex spatial distribution, 
we perform a spatial average for all the positions with 
crb,f > 1 and also perform a time average from t = 90 to 
t = 120 to get the dby ~ 1.16. As a result, we derive 
77 ~ 35.7% based on the analytical calculation (Eq.([2Q|)). 
This is well consistent with the Eem dissipation efficiency 
calculated directly from the energy evolution of the sim¬ 
ulations using Eq.(13), as shown in Figure 


3.6. Summary for this section 

In this section, we revealed a collision-induced strong 
reconnection process with the EMF energy dissipation ef¬ 
ficiency about 35%, which is resolution insensitive. The 

® However, the misalignment Xg itself does have a direct effect 
on the Eeui dissipation efficiency due to the different field config¬ 
urations around the initial contact surface. See details in Section 
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Fig. 7.— The contour cuts corresponding to the maximum y-component of the outflow velocity (Vy) in the TZ-plane (x=0) for different 
resolutions. The maximum values of Vy are about 0.45c, 0.55c, 0.75c for the resolutions of 256^, 512^, and 1024^, respectively. The aspect 
ratio becomes smaller for a higher resolution. 



Fig. 8.— The evolution of the held lines during the collision process. The two blobs merge into one larger blob, forming a “oo”-shaped 
held line conflguration at the flnal quasi-steady stage of the evolution. 


outflow can locally become super-Alfvenic during the 
initial strong compression stage. The outflow velocity 
can potentially become relativistic in higher resolution 
simulations and generate multi-orientation mini-jets in a 
global PFD jet. 

4. EXTENDED PARAMETER SPACE STUDIES 

In Section ^ we find significant EME energy dissipa¬ 
tion (about 35%) facilitated by collision-driven magnetic 
reconnection. Based on the above analyses, we expect 
that some parameters may affect the results. Eirst, the 


(75 evolution controls the dissipation efficiency. The 
simulation results may then depend on the initial value 
(7b,i. Second, the initial misalignment Xg gives different 
magnetic field configurations around the contact surface 
which may control the fraction of the free energy that 
can be released due to the reconnection processes. Next, 
different initial relative speed (kinetic energy) between 
the two blobs define the strengths of the initial collision- 
driven effect, so that it may be another factor to effect 
the conclusion. In addition, the initial displacement 
controls the delay of the collision. It is also interesting 
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to investigate whether results depend on this parameter. 
Furthermore, the blobs undergo a significant expansion 
during the early “self adjustment” phase before estab¬ 
lishing a balance between the magnetic pressure force 
and the gas pressure force. Different background pres¬ 
sures and densities are therefore interesting input pa¬ 
rameters that may affect the results. Finally, it is un¬ 
known whether the toroidal-to-poloidal ratio parameter, 
a, plays a role to define the dissipation efficiency. 

In this section, we perform a series of extended parame¬ 
ter studies to investigate the role of above-mentioned var¬ 
ious parameters in defining the F^em dissipation efficiency 
and reconnection outfl ow properties. From the resolution 
study in Section §3.41 we And that the F^em dissipation 
efficiency is similar when the resolution is ^ 256^. Since 
we are exploring a large parameter space, in order to 
reduce the simulation time, we use the 256^ resolution 
in all the simulations presented in this section. The gen¬ 
eral approach is that we only modify one parameter from 
the example simulation in each subsection, in order to ex¬ 
plore the effects of that parameter. Below we explore the 
effect of following parameters in turn: the initial a value 
inside the blobs (crb,i), the initial misalignment between 
the center of two blobs in X direction (x^), the initial ve¬ 
locity (Vb,z) or relative Lorentz factor (Frei) of two blobs 
in Z direction, the initial distance between the center 
of two blobs in Z direction (zd), the uniform gas pres¬ 
sure (P), the background density (pbkg), the toroidal-to- 
poloidal magnetic held ratio a, and the adiabatic index 
(7)- 


4.1. Initial a\)^i of the blobs 


The crb,i parameter affects the evolution and the ini¬ 
tial expansion of the blobs, so in this subsection we study 
the effect on the Pem energy evolution and dissipation 
efficiency for different <Jb,i values. We set a higher crb,i 
value by increasing the normalization parameter Pb,o- 
The parameters we used are listed in Table Here the 
parameters besides crb,i and Pb,o are the same as the 
example model. We choose three different crb,i values. 
Here, due to the more significant expansion with a in¬ 
creasing crb,i, we enlarge our simulation box to 30^ and 


also increase the resolution to 384^ to ensure the same 
absolute spatial resolution as the following sub-sections. 
The upper panel of Figure shows the Pem evolution of 
these three models. Similar to the example model, we 
also show the non-collision case for each model to calcu¬ 
late the Pem dissipation efficiency. We And that with a 
higher crb,i value, the initial steep decay of Pem caused 
by the expansion is more significant due to the stronger 
outward magnetic pressure force. The Pem dissipation 
efficiency in the Anal quasi-steady phase, on the other 
hand, is rather similar for diflerent ab.i values (lower 
panel of Figure 1^. In order to understand this result, we 
also calculate the dby values in the quasi-steady phase, 
and use Eq.(20) to cal culate 77 independently (see the 
method in Section §3.5| ). The results are shown in Table 
m The calculated efficiencies have no obvious diflerence 
among diflerent cases, since a larger crb,i corresponds to a 
slightly larger so that 77 does not change signiAcantly 
according to Eq.(20). This result is consistent with the 
efflciency calculatea directly from the energy evolution 
of the simulations in the lower panel of Eigurej^ 

One interesting result is that crb,i and crb,f values show 


TABLE 3 

The (7b,i-VARYING MODELS 


Model name 

*^b,i 

^b,0 

a 

l^b,z| 

P 

Pbkg 

Zd 

Xs 

Model:(7b,i8 

8 

VTtt 

3 

0.3c 

(M 

1 

o 

1—1 

O 

1 

4.4 

1.0 

Model:c7b,il6 

16 

\/Xk 

3 

0.3c 

(M 

1 

o 

1—1 

o 

1 

4.4 

1.0 

Model:(7b,i24 

24 


3 

0.3c 

10-2 

10-1 

4.4 

1.0 
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Time 


Fig. 9. — The -Eem evolution of three crb,i models: (7b,i = 8 (red 
group), (7b,i = 16 (green group), and C7b,i = 24 (blue group). As the 
example model, we also show the non-collision cases corresponding 
to each of them to calculate the Eem dissipation efficiency. The 
Bern dissipation efficiency in the final quasi-steady phase is similar. 


TABLE 4 

f^b,f - f^b,i RELATION AND THE ANALYTIGAL VS. NUMERIGAL 
EFFIGIENGIES. 


^b,i 

f^b,f 

Efficiency 

(analytical) 

Efficiency 

(numerical) 

8 

1.16 

35.7% 

33.3% 

16 

1.33 

37.3% 

34.4% 

24 

1.49 

36.4% 

34.7% 


a good linear relationship (EigurefTol), which can be Atted 
by 


crb,f = 0.02crb,i + 1.0. (22) 


Physically, this equation does not apply for crb,i < 1. 
Right now the range of the crb,i is relatively small limited 
by the code capability. It is valuable to perform a more 
detailed study for a larger range of crb,i in the future to 
test this interesting and important relationship and to 
study the physical mechanism of this behavior. 


4.2. Initial misalignment between two blobs (xg) 

The initial misalignment affects the magnetic held con- 
Aguration around the contact surface during the collision 
and also the rotation property, so in this subsection we 
study the effect of Xg on the F^em energy evolution and 
dissipation efflciency. The parameters we used are listed 
in Table Besides Xg, other parameters are the same 
as th e example model. From the analyses in Section 
§3.5[ the typical radius of one blob after expansion is 
r « 3ro = 3.0. In order to make a relatively signiAcant 
collision, Xg should be smaller than 2r ~ 6. We choose 
seven diflerent values of Xg in the simulations (Table 
Among these models, “Modehxgl.O” is the same as the 
example model with resolution 256^. 

The upper panel of Figure pT] shows the Egm evolution 
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1.5 

X ^b,'f 

- Linear fitting: Ob f=0.02Gb j+1.0 


1.4 



1.3 






1.2 



1.1 




5 10 15 20 25 30 

^b.i 

Fig. 10.— The relationship between (7b,f and (7b,i. The red cross 
points are the calculated results of crb,f corresponding to different 
(7b,i from our simulations. The black dashed line is the linear fitting 
result. 

of these seven models. We find that, with increasing 
the E’en! dissipation efficiency first increases, and then 
decreases after reaching the maximum efficiency around 
Xg ~ 3.0. Such a behavior can be understood using the 
two lower panels in Figure The 2D cuts of By 
show the directions and strength distributions of these 
two components of the magnetic field. Due to the ini¬ 
tial expansion, the actual size of these configuration dur¬ 
ing collision would become about three times larger than 
the initial cuts. We consider the anti-parallel regions 
as the free energy source for reconnection-induced dis¬ 
sipation. When Xg = 0, only the B^ component can 
be reconnected. With an increasing Xg from 0 to 3.0, 
the anti-parallel region of B^ becomes smaller while the 
anti-parallel region of By increases. Since the strength 
of By is larger than B^^ the total dissipated magnetic 
energy becomes larger, which increases the dissipation 
efficiency to a higher value. When Xg = 3.0, the anti¬ 
parallel region of By reaches the maximum so that the 
maximum efficiency is achieved. After this critical point, 
the anti-parallel region of By decreases with increasing 
Xg^ which leads to a decrease in efficiency. Finally, when 
Xg becomes larger than the size of the two blobs (6.0), 
the F^em evolution is nearly the same as the non-collision 
case due to the lack of significant collision between the 
two blobs. 

Another important feature is that the change of effi¬ 
ciency as a function of Xg is not linear. In fact, in our 
studied the cases only the two extreme cases = 0 
and Xg = 7.0 have significantly lower efficiencies com¬ 
pared with the other cases. Even if there is a very small 
misalignment, e.g. Xg = 0.002, the efficiency could be¬ 
come significantly larger than the Xg = ^ case. Inspecting 
the energy evolution plot (Fig|TT|, one can see that the 
main difference comes from the~~^‘normal decay” phase. 
A smaller Xg would delay the “normal decay” phase to 
a later time, whereas there is no “normal decay” phase 
if Xg = 0. In order to understand this feature, we draw 
Figure The three panels in the first row are the 2D 
contour cuts of B‘^ jAn in the XZ-plane (y=0), corre¬ 
sponding to t = 0, t = 80, and t = 120 for “ModehxgO.O”, 
respectively. The following two panels in the second row 
are the 2D contour cuts of the ^-component of the out¬ 
flow velocity {Vy) in the TZ-plane (x=0), correspond¬ 
ing to t = 80 and t = 120 for “ModehxgO.O”. The 


following five panels are the corresponding cuts for the 
“Model:Xs0.002”. One can see that for the complete 
alignment case (“ModehxgO.O”), the two blobs merge to¬ 
gether and enter the quasi-steady phase without an ob¬ 
vious reconnection feature. For “Model:Xs0.002”, on the 
other hand, the merged blobs start to rotate at t ^ 80, 
which just corresponds to the starting time in the energy 
evolution plot (Fig pT| when the two models become dis¬ 
tinctly different. The last panel of Vy in Figureshows 
the feature of the reconnection-driven outflows during 
the “normal decay” phase of “Model:Xs0.002”, which 
is not seen in “ModehxgO.cQ From these analyses, we 
can draw conclusion that, besides the initial anti-parallel 
region caused by Xg^ the collision-triggered rotation is 
another important process to change the magnetic field 
configuration and dissipate more magnetic energy. This 
rotation-driven dissipation seems to correspond to the 
“normal decay” phase in the energy evolution plot. 


TABLE 5 

The Xs-varying models 


Model name 

*^b,i 

^b,0 

a 

Ab.zl 

P 

Pbkg 

Zd 

Xg 

ModeharsO.O 

8 

VVt 

3 

0.3c 

10-2 

10-1 

4.4 

0 

Model:a:s0.002 

8 

VTtt 

3 

0.3c 

10-2 

10-1 

4.4 

0.002 

Model:a:s0.02 

8 

\/47r 

3 

0.3c 

(M 

1 

o 

1—1 

10-1 

4.4 

0.1 

Modeha^sO.l 

8 

VVk 

3 

0.3c 

(M 

1 

o 

1—1 

10-1 

4.4 

1.0 

Modeha^sl.O 

8 

\/47r 

3 

0.3c 

(M 

1 

o 

1—1 

10-1 

4.4 

3.0 

Modehxsl.h 

8 

VTr 

3 

0.3c 

(M 

1 

o 

1—1 

10-1 

4.4 

5.0 

Modehxsl.h 

8 

VTr 

3 

0.3c 

(M 

1 

o 

1—1 

10-1 

4.4 

7.0 
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Fig. 11.— The Eem evolution of six Xg models: Xg = 0 (red), 
Xg = 0.002 (green), Xg =0.1 (blue), Xg = 1.0 (magenta), Xg = 3.0 
(cyan), Xg = 5.0 (yellow), and Xg = 7.0 (gray). Similar to the 
example model, we also show the non-collision cases for each model 
to calculate the Eem dissipation efficiency. Since the non-collision 
case does not depend on Xg , there is only one non-collision evolution 
case (dashed line). Besides the Xg7.^ model which no collision is 
observed. Only the “xsO.O” model has obvious difference in Eem 
dissipation efficiency, which is only about half of the other models 
in the final quasi-steady phase. 


4.3. Initial relative Lorentz faetor (Treiy^ between the 
two blobs 

A larger initial velocity Vb,z means a larger initial ki¬ 
netic energy of the two blobs, which would provide a 

^ The outfiows look like asymmetric on that 2D cut, but the 3D 
configuration is more complex. 
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Fig. 12.— Selected cuts for model “Modelia^sO.O” (first two rows) and “Modelia^s0.002” (last two rows). For each model, we draw 
the 2D contour cuts of jAir in the XZ-plane (y=0) at t = 0, t = 80 and t = 120 (higher row) and the 2D contour cuts of the y- 
component of outflow velocity (Vy) in the yZ-plane (x=0) at t = 80 and t = 120 (lower row). The results show that even if there is a very 
small misalignment, the merged blobs would start to rotate at some critical point and trigger additional reconnection-facilitated magnetic 
dissipation. 


larger driving force initially and also effectively decrease factor can be calculated as 
crb,i. In this subsection, we study the effect of Vb,z in de¬ 
tail. Since we are testing Vb,z in the relativistic regime, 
we adopt the relative Lorentz factor (Frei) between two 

blobs as the varying parameter. The relative Lorentz Trei = 2r^ ^ ~ 


(23) 
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where rb,z = (1 ~ z/^^) Lorentz factor of 

each blob. 

The parameters we used are listed in Table The 
parameters except Trei and crb,i are the same as the ex¬ 
ample model. The effective change of crb,i is a conse¬ 
quence of changing Trei, since we keep B as observed in 
the center-of-mass frame as constant. We choose three 
different Trei values. With an increasing Trei, the frac¬ 
tion of the initial kinetic energy becomes larger, so that 
crb,i effectively decreases. The upper panel of Figur e p!3| 
shows the £^em evolution of these three models. As i Vei 
increases, due to the stronger dynamic process, the ini¬ 
tial “self adjustment” phase becomes more significant. 
In particular, for Trei = 18.8, the fast motion generates 
a strong shock into the surrounding medium. The shock 
compresses the magnetic field even further, so that £^em 
reaches an even higher value initially. In addition, from 
the lower panel of Figure we find that the F’em dissi¬ 
pation efficiency in the final quasi-steady phase increases 
with increasing Frei- This is because a higher initial Frei 
carries a larger kinetic energy and gives a stronger initial 
collision-induced compression and reconnection-induced 
dissipation. 


TABLE 6 

The Prel-VARYING MODELS 


Model name 

c^b,i 

-6b,0 

a 

Trel 

P 

Pbkg 

Zd 

Xs 

Model:F^eH-2 

9 

VTtt 

3 

1.2 

(M 

1 

o 

1—1 

O 

1 

4.4 

1.0 

Model:F^ez5.6 

6 

y/pK 

3 

5.6 

10-2 

10-1 

4.4 

1.0 

Model:F^eP8.8 

3 

y/pK 

3 

18.8 

10-2 

10-1 

4.4 

1.0 
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Fig. 13. — The Eem evolution of four Frei models: Frei = 1-2 
(red), Frei = 5.6 (green), and Frei = 18.8 (blue). Similar to the 
example model, for each model, we also show the non-collision 
case (dashed) to calculate the Eem dissipation efficiency, the Eem 
dissipation efficiency is larger for the model with a larger Frei- 

4.4. Initial distance between two blobs (zd) 

The initial distance between the center of two blobs 
in z direction (zd) controls the delay of the collision. In 
our example case, Zd is relatively small and the collision 
happens around the middle stage of the initial expan¬ 
sion. It is valuable to study the effect on the energy 
evolution and dissipation efficiency when we set a larger 
Zd to delay the collision time to a later stage. The pa¬ 
rameters we used are listed in Table Here the param¬ 
eters except Zd are the same as the example model. We 


choose three different Zd values. The “Model:2;c^4.4” uses 
Zd = 4.4 which is the same as the example model with 
resolution 256^, and the other two models, “Modeh^^^^G” 
and “Model:Zc^8”, have Zd = 6,8, respectively. The up¬ 
per panel of Figure [T^ shows the F^em evolution of these 
three models. The collision times of these three different 
Zd models are designated to be around the middle and 
late stages of initial expansion as well as after the initial 
expansion, respectively. We find that the reconnection- 
triggered F^em dissipation process is systematically de¬ 
layed when Zd becomes larger. However, the efficiency at 
the final quasi-steady phase reaches a similar value (lower 
panel of Figure jT3 in all three models. This result sug¬ 
gests that the dissipation efficiency does not sensitively 
depend on the stage of blob evolution. 


TABLE 7 

The 2;c^-varying models 


Model name 

f^b,i 

6b,0 

a 


P 

Pbkg 

Zd 

Xs 

Mode\:zd4:A 

8 

y/pK 

3 

0.3c 

o 

1 

to 

10-1 

4.4 

1.0 

Model:2;c^6 

8 

y/pK 

3 

0.3c 

(M 

1 

o 

1—1 

10-1 

6 

1.0 

Model:2;c^8 

8 

pPPr 

3 

0.3c 

(M 

1 

o 

1—1 

10-1 

8 

1.0 
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Fig. 14.— The Eem evolution of three Zd models: Zd = 4.4 (red), 
= 5 (green), and Zd = S (blue), which correspond to the collision 
happening at different stage of the blob evolution. Similar to the 
example model, the non-collision case (black dashed, the same for 
all three models) is plotted for comparison. Although the reconnec¬ 
tion facilitated dissipation is delayed when Zd increases, the energy 
level and the Eem dissipation efficiency in the final quasi-steady 
phase are similar. 


4.5. Background gas pressure (P) 

The blobs undergo a significant initial expansion to 
establish a balance between the magnetic pressure and 
the ambient gas pressure. Different background pressure 
(P) and density (pbkg) would affect these processes. In 
this and next subsections, we study the effect of P and 
Phkg on the Pern energy evolution and dissipation effi¬ 
ciency. The parameters for the P-varying models are 
listed in Table Here the parameters other than P 
are the same as the example model. We choose three 
pressure values. The “Model:P-2” is the same as the 
example model with resolution 256^, and the other two 
models have a lower (Model:P-I) or higher (ModehP- 
3) pressure. The upper panel of Figure shows the 
Pem evolution of these three models. The initial “self 
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adjustment” phase due to the expansion is more signifi¬ 
cant when pressure becomes lower, since for a lower back¬ 
ground pressure it takes longer for the magnetic blob to 
expand before reaching a balance with the ambient gas, 
and vice versa. Even with very different blob dynamics 
for different pressure values, the difference of dissi¬ 
pation efficiency in the final quasi-steady phase is not so 
large as the evolution itself (lower panel of Figure 


15), which means that the efficiency is relatively insen¬ 


sitive to the expansion process. This is probably due to 
the fact that the initial expansion phase with different 
background pressure values affects the evolution of both 
the collision and non-collision cases to similar degrees, 
so that the relative difference (efficiency) does not have 
a significant change. 


TABLE 8 

The P-varying models 


Model name 

f^b,i 

^b,0 

a 

IVz| 

P 

Pbkg 

Zd 

Xs 

Model:P-l 

8 

VTr 

3 

0.3c 

10-1 

10-1 

4.4 

1.0 

Model:P-2 

8 

VTr 

3 

0.3c 

10-2 

10-1 

4.4 

1.0 

Model:P-3 

8 

VTr 

3 

0.3c 

10-3 

10-1 

4.4 

1.0 
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Fig. 15. — The Pem evolution of three pressure models: P = 
10“^ (red), P = 10“^ (green) and P = 10“^ (blue). Similar to 
the example model, we also show the non-collision cases (dashed 
lines) corresponding to each pressure model to calculate the Pem 
dissipation efficiency. Although the initial “self adjustment” phase 
due to the expansion becomes more significant when the pressure 
goes to lower values, the difference of efficiency in the final quasi¬ 
steady phase is relative small. 


4.6. Background density (phkg) 

The parameters for pbkg-varying models are listed in 
Table Here most of the parameters are similar to the 
example model. However, in order to make the code sta¬ 
ble when applying a smaller background density, we have 
to increase the uniform gas pressure value to 10“^. We 
choose two density values. The upper panel of Figure p!6| 
shows the F^em evolution of these two models. There is 
no significant difference between these two models. This 
is understandable. Since the force balance is mainly con¬ 
trolled by the background gas pressure rather than den¬ 
sity, varying the background density does not lead to 
significant change in the result. 


TABLE 9 

The Pbkg-VARYING MODELS 


Model name 

*^b,i 

^b,0 

a 

l^b,z| 

P 

Pbkg 

Zd 

Xs 

Model:pbkg-l 

8 


3 

0.4c 

o 

1 

O 

1 

4.4 

1.0 

Model:pbkg-3 

8 


3 

0.4c 

o 

1 

10-3 

4.4 

1.0 
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Fig. 16. — The Pem evolution of two density models: Pbkg = 
10“^ (red), and Pbkg = 10“^ (green). Similar to the example 
model, the non-collision cases are also plotted (dashed). There is 
essentially no difference between the two models. 


The a parameter (introduced in Eq.(|^) defines the 
ratio between the toroidal and poloidal components of 
the initial magnetic field configuration. In the example 
model, we adopted a = 3, which means that the toroidal 
flux roughly equals to the poloidal flux. On the other 
hand, the central engine of GRBs (and probably AGNs 
as well) is likely rapidly rotating, so that the toroidal flux 
may be (much) larger than the poloidal flux and a > 3. 
In this subsection we study the effect of a on the Fem en¬ 
ergy evolution and dissipation efficiency. The parameters 
we used are listed in Table pT} Here the parameters other 
than a are the same as the example model. We choose 


two a values (3 and 8). The upper panel of Figure 17 
shows the Fem evolution of these two models. The init^ 
“self adjustment” phase due to the expansion is more sig¬ 
nificant when a becomes higher, since th e net magneti c 
pressure becomes larger when a increases ( Li et al.|2QQ6 ). 
This requires more expansion of the blobs before reach- 
ing the pressure balance with the ambient gas. The lower 
panel of Figure shows the F^em dissipation efficiency of 
these two models. The model with a higher a value has 
a slightly higher efficiency. This is probably because the 
reconnections mainly come from the toroidal component. 
If a is much greater than 3 in realistic astrophysical sys¬ 
tems (e.g. in GRBs and AGNs), the efficiency derived 
here can be regarded as a safe lower limit. 


TABLE 10 

The a -varying models 


Model name 


^b,0 

a 

l^b,z| 

P 

Pbkg 

Zd 

Xs 

ModebaO 

8 

\/47r 

3 

0.3c 

10-2 

10-1 

4.4 

1.0 

Model:Q:8 

8 

vTtt 

8 

0.3c 

10-2 

10-1 

4.4 

1.0 


4.7. a value 


4.8. Adiabatic index 
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Fig. 17.— The £^em evolution of two a models: a = 3 (red) and 
o; = 8 (green). Similar to the example model, the non-collision 
cases for each model (dashed) are also plotted. The -Fem dissipation 
efficiency is slightly larger for the model with a larger a. 


Finally, in all above simulations we have used a sim¬ 
ple uniform adiabatic index 7 = 5/3, since most of the 
regimes are at most mildly relativistic. Nonetheless, in 
some high temperature regimes in the reconnection out¬ 
flows, the adiabatic index may be close to the relativistic 
limit 7 = 4/3. In principle, for a more accurate study, we 
need to calculate the adiabatic index between 5/3 to 4/3 
based on the temperature distribution and time evolution 
cell by cell. In this subsection, based on the simplified 
uniform adiabatic index model, we compare the differ¬ 
ence between these two limits: 7 = 5/3 and 4/3. We 
name them as “Modehadi 5/3” and “Modehadi 4/3”, 
respectively. The parameters that we used are listed 
in Table [ll| The “Modehadi 5/3” is just the exampL 
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case with resolution 256^. The upper panel of Figure 
shows the F^em evolution of these two models. Similar to 
the example model, we also show the non-collision case 
for “Modehadi 4/3”. We find that there is only a slight 
difference between these two models. We therefore con¬ 
clude that the simple uniform adiabatic index model with 
7 = 5/3 catches the essence of the collision and magnetic 
dissipation physics discussed in this paper. 



0 20 40 60 80 100 120 

Time 


Fig. 18.— The Eem evolution of two 7 models: 7 = 5/3 (red) and 
7 = 4/3 (green). The non-collision cases (dashed) for both models 
are also shown. The -Fem evolution and dissipation efficiency are 
similar for these two models. 


4.9. Summary for this section 

In this section, we have done a series of extended pa¬ 
rameter studies. We find that the F^em dissipation effi¬ 


ciency is relatively insensitive to the variations of most 
parameters compared with the F^em evolution itself. We 
conclude that the two conclusions drawn in Section ^ 
namely, a high collision-triggered magetic dissipation ^ 
ficiency and the existence of reconnection-facilitated out¬ 
flow minijets, are robust. 


5. MULTIPLE COLLISIONS 

So far we have only studied collisions between two high¬ 
er blobs. In reality, multiple collisions between several 
or even tens of blobs may occur in GRB/AGN jets, so 
that a much more complex configuration in the collision 
regions with multiple reconnection layers and outflows 
can be achieved (a s envisaged in the IGMART model, 
Zhang fc Yan||2Qll ). Simulating multiple collisions with 
random blobs are technically heavy duty and require ex¬ 
tended simulation efforts. Nonetheless, in this section, 
we present a preliminary four-blob interaction simula¬ 
tion as the first step towards a more realistic/complicated 
study of multiple-blob collisions. 

The size of the simulation box is 20 x 20 x 40 with res¬ 
olution 256 X 256 X 512, which means that the grid size is 
the same as the example case of two blobs with resolution 
256^ in Section ^ The parameters used in the simula¬ 
tion are listed in Table Here most of the parameters 
are the same as the example model. The different ones 
include the initial velocities of the four blobs and the dis¬ 
tances between blobs. For the four blobs along the -\-Z 
direction (bulk motion direction of the global jet) in the 
simulation frame (center-of-mass frame of the blobs), the 
velocities are 0.8c, 0.3c, —0.3c, and —0.8c, respectively. 
The initial distances between them are 5, 15, and 5, re¬ 
spectively, which means that the inner two blobs would 
collide with their nearby outer neighbors first to form 
two merged knots, before the two knots collide again. 
Figureshows the 3D contour plots of the current, and 
the corresponding 2D contour cuts of the ^-component 
of the outflow velocity(U^) in the YZ plane. At t = 28 
the first collisions between the two pairs of outside blobs 
form two strong reconnection regions which are similar 
to the example model in ^ At later times, these two 
post-collision knots collide again and form a third strong 
reconnection region in the middle, and the original two 
reconnection regions also continuously evolve with time. 
These three reconnection regions form a more complex 
configuration than the collision with only two blobs as 
studied in ^ 

The upper panel of Figure shows the Eem evolu¬ 
tion of this case. For comparison, we also show the 
non-collision case and the two-blob collision case (with 
same resolution), and calculate the Eem dissipation ef¬ 
ficiency. We find that the Eem dissipation efficiency in 
the final quas i-steady phase is around 40% (lower panel 
of Figure ]2Q|), which is higher than the example two- 
blob case with resolution 256^ (lower panel of Figure]^. 
This suggests that multiple collisions can facilitate mr- 
ther reconnection-triggered magnetic dissipations, mak¬ 
ing the system reaching a higher Eem dissipation effi¬ 
ciency. 


6 . GONGLUSIONS AND ASTROPHYSIGAL 
APPLIGATIONS 

In this paper, using a 3D SRMHD code, we carried out 
a series of simulations to study collisions between high- 
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TABLE 11 

Two 7 MODELS 


Model name 

7 


^b,0 

a 

l^b,z| 

P 

Pbkg 

Zd 

Xs 

Modekadi 5/3 

5/3 

8 


3 

0.3c 

10-2 

10-1 

4.4 

1.0 

Modekadi 4/3 

4/3 

8 

VTr 

3 

0.3c 

10-2 

10-1 

4.4 

1.0 


TABLE 12 

The four-blob collision model 


Model name f7b,i ^b,o 

a Eb,z P 

Pbkg 

Zd Xs 

Model: 4 blobs 8 VTtt 

3 0.8c/0.3c/-0.3c/-0.8c 10-2 

10-1 

5/15/5 1.0 



15 -10 -5 0 5 10 15 


Eig. 19.— Some results of the four-blob collision simulation. The three left panels are the 3D contour plots of the currents at t = 28, 
t = and t = 64, respectively. The two right panels are the corresponding 2D contour cuts of Vy in the YZ plane (x=0) at t = 28 and 
t = 64. The black lines indicate the corresponding positions of strong current layers and outflows. At least three strong reconnection layers 
with different directions are formed during the collision process. 


a magnetic blobs. Through a detailed example simula¬ 
tion and an extended parameter space survey, we have 
reached the following robust results: 


parallel regions in the contact surface of the two 
blobs. As long as a small offset exists, significant 
dissipation is facilitated. 


• Collisions trigger significant EMF energy dissipa¬ 
tion. Detailed analyses of the numerical data dur¬ 
ing different stages of the collision process sug¬ 
gest that such dissipation is facilitated by collision- 
induced magnetic reconnection. The efficiency of 
Fern dissipation in our simulations is around 35%, 
which is insensitive to the numerical resolution and 
several initial condition parameters, such as crb,i, 
Trei, Zd, P, and Pbkg- It depends on the impact 
parameter Xg, which defines the area of the anti- 


• Our simulations suggest that the collision process 
is essentially inelastic. Even though there is some 
kind of bouncing back in the early stage of the col¬ 
lision evolution, the strong reconnection effect in 
the contact surface efficiently dissipates the mag¬ 
netic energy and reduces the magnetic pressure. 
As a result, the two high-cr blobs merge into one 
larger blob with a “oc”-shaped m agne tic configu¬ 
ration (see more details in Section ^3.5). Assuming 
complete inelastic collision, an estimated dissipa- 
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Fig. 20.— The Bern evolution of four-blob model compared with 
the two-blob model at resolution 256^. For comparison, the non- 
collision model is also shown. The -E’en! dissipation efficiency in 
the final quasi-steady phase for the four-blob model is significantly 
higher than that of the two-blob model. 


tion efficiency (Zhang & Yan[2011) is found consis¬ 
tent with the efficiency derived from the numerical 
data. 


• From our simulations, we find that magnetic re¬ 
connection events can induce relativistic, multi¬ 
orientation outflows. Even in two-blob simulations, 
as long as an offset exists {xs 7 ^ 0), 3D outflows are 
formed. For four-blob collisions, significant out¬ 
flows exist in three distinct regions. These out¬ 
flows would mimic “mini-jets” as invoked in the 
astrophysical models of GRBs and AGNs. The 
maximum outflow velocity (V^ut) in owe simula¬ 
tions is only 0.75c. However, due to the signifi¬ 
cant resolu tion- dependence behavior as described 
in Section §3.4[ we still have not reached the con¬ 
vergence for the outflow velocity, so that that value 
is only the lower limit of hbut- In the simulations, 
we found that Tp^ t can reach and even exceed local 
Ta and T^s (^3.3), both are relativistic numbers if 
cr > 1 . Also a larger Trei tends to give a larger Vout 
value. In principle, with a high-resolution simula¬ 
tion and for a large crb,i and Trei, an even larger 
mini-jet Lorentz factor is achievable. 


• We found an interesting linear relationship between 
the initial (crb,i) and the flnial (ab.f) values of the 
a parameter of the blob (Eq.([2^). The range of 
crb,i we have explored is not very large due to the 
code capability constraint. It is valuable to study 
this intriguing behavior in a larger range of crb,i in 
the future. 


• Our preliminary simulations of multiple collisions 
among multiple high-a blobs suggest that the 
collisions would give rise to more complex con¬ 
figurations of the reconnection layers and multi¬ 
orientation outflows, with a higher EMF energy 
dissipation efficiency. This suggests that the multi¬ 
ple collisions of many high-cr blobs can potentially 
generate many mini-jets with relatively random di¬ 
rections, as required by some theoretical models of 
astrophysical jets. 


These numerical simulations have profound implica¬ 
tions to understand astrophysical jets, such as GRBs, 


AGNs, X-ray binaries. Grab nebula, and so on. In the 
following, we discuss their direct applications to GRB 
and AGN models. 


6.1. GRBs 


As w e ment ioned in the introduction section, |Zhang| 
& Yan (2011) proposed the ICMART model to inter- 
pret the prompt emission of GRBs. This model in¬ 
vokes collision-induced magnetic dissipation of moder¬ 
ately high-cr blobs, which is the motivation of our sim¬ 
ulations. The ICMART model was suggested to have 
several salient features that can potentially interpret var¬ 
ious observations not easy to interpret within the MDF 
internal shock models. Our simulations verified several 
as sumptions/speculatio ns adopted in the original model 
of iZhang fc %n\ (|2QTT]L 



jZhang Yan (2011) claimed that ICMART pro¬ 
cesses snould have a significantly higher energy dissipa¬ 
tion efficiency than internal shock s, which is more con¬ 


sistent with the GRB obs ervations (Panaitescu & Kumar 
IZhang et al.||2007|). They assumed that once IC 

T 


2002 

MAIIT' is triggered', the two colliding shells would merge 
completely in an inelastic collision. The a values of the 
two shells/blobs drop significantly from an initial value 
to a much lower final value. Energy and momemtum 
conservations suggest that the energy dissipation effi¬ 
ciency is high, up to 10 s of percent, depending on the 
final cr value of the merged blob. If crb,f ^ 1, they found 
that the efficiency is close to 50%. Our detailed simula¬ 
tions verified all these assumptions/speculations. Indeed 
significant magnetic dissipation occurs due to collision- 
induced magnetic reconnection. The collision process is 
essentially inelastic, and the energy dissipation efficiency 
is indeed high, which is ^ 35% in for two-blob collisions 
and ^ 40% for four-blob collisions. One surprising result 
is that the final value dby is linearly correlated with the 
intial value crb,i (Eq.([2^), so that the efficiency does not 
sensitively depend on (Jby. More studies are needed to 
reveal the underlying physics of this correlation. 

Second, the IGMART model invokes the central engine 
activities to interpret the broad pulses in the GRB light 
curves, but requires the existence of mi ni-jets to account 
for th e rapid variability component. | Zhang fc Zhang] 
(2014) used this concept to perform a series of Monte 
Carlo simulations and reproduced a range of highly vari¬ 
able light curves with both s low and fast com ponents as 
seen in observational data (|Gao et al.|| 2012 |). The re- 
quire d Lorentz factor of th e mmi-jets is in the range of 
2-15 (Zhang & Zhang|2014). In our simulations, the out¬ 
flows of reconnection layers can reach mildly relativistic 
speed. From the orientation point of view, one major 
reconnection current layer between two colliding blobs 
already generates multi-orientation outflows (see Figure 
[^, in addition to systematical global rotation and twist 
one to the slightly initial misalignment. Furthermore, by 
invoking four-blob collisions, we clearly find three ma¬ 
jor reconnection layers with different directions. Each of 
them has their own 3D outflow systems similar to two- 
blob collision cases, which gives a more complex space- 
time distribution of the outflow directions. In realistic 
systems, collisions of tens of blobs would lead to more 
complicated 3D mini-jet structure, which would account 
for observed GRB light curves. It is possible in much 
smaller scales not resolved by the current simulations. 
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perturbations may induce turbulent reconnections, which 
may make even smaller mini-jets in the current outflows. 
Dedicated local simulations are needed to verify or refute 
such a speculation. With the current global simulations, 
one is confident that even without turbulence, collision- 
induced reconnection layers can already generate large- 


local magnetic reconnections in a global, Poynting-flux- 
dominated jet. These local reconnections generate lo¬ 
cal outflows or mini-jets with a comoving Lorentz fac¬ 
tor around a few. Our simulations give an alternative 
process to trigger the local reconnections by consider¬ 
ing ICMART events, i.e. collisions among magnetic 


scale mini-jets in the bulk jet of a GRB, which would 

knots/blobs inside the global jet. 

Since the knots in 

give interesting variability features in the light curves. 

AGNs have already been observed ( 

Marscher et al.|2002 

Einally, our simulations show significant evolution of 

Chatterjee et al. 

2009 Doi et al. 

201l|), the colnsions 


event. From Figure we can see that during the early 
collision-driven reconnection phase (e.g. the “plateau” 
phase), the strengths of and By components de¬ 
crease and the Bz component increases, which signifi¬ 
cantly changes the magnetic field configuration while still 
keeping a relative ordered magnetic configuration in a 
relatively short time duration. The behavior may poten¬ 
tially interpret the significant change of the polarization 
angle duri ng the prompt emissio n phase observed in GRB 
100826A ( Yonetoku et al.|[2Qll ). 


6.2. AGNs 

Some blazars show very fast TeV flares whose dura¬ 


tions are only sev eral minutes (Aharonian et al. 2007 
Albert et al.||2007 ). This duration is much shorter than 


the light crossing time for the entire system, which means 
that emission comes from a small local region. The re¬ 
quirement of emitting TeV photons also demand a much 
lar ger Lorentz factor in the emission region (greater than 
50, Begelman et al.||2008 Mastichiadis fc Moraitis||2008 ) 
than what is i nferred tor the bulk motion (typically 
small er than 10, Giroletti et al.|[^04| [Finer fc Edwards 
2004). A successful model to interpret t he observations 
i s the “jets in a jet” model proposed by [Giannios et al. 
(2009). This model invokes current-instability triggered 


would very likely happen, which trigger the local recon¬ 
nections and generate the mini-jets as needed in their 
model. 

We can also roughly estimate the time scale using our 
simulation resu l ts an d the parameters in the model of 


Giannios et al. (2009). They estimated that the typical 
size of the blob in the rest frame of the blob is around 
lO^^cm. Since the Lorentz factor of the blob in the co¬ 
moving frame of the global jet is equal to 10 (assuming 
(j = 100), the size of the global jet in the comoving frame 
is about lO^^cm, which can be treated as Lq in our Ta¬ 
ble [H Thus one time unit in our simulations can be 
normalized as to = Lq/c ^ 3 x lO^s. The duration of 
the reconnection-facilitated energy dissipation is about 
50 time units in our example case (see Figure from 
t=30 to t=80), which can be translated to about lO^s 
in the rest frame of the global jet. The duration in the 
observer frame is ^ lO^s/Fj ^ lO^s/10 = lO^s, which is 
very close to the observed durations of the flares. 
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